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APPROXIMATE  DIRICHLET  BOUNDARY  CONDITIONS  IN 
THE  GENERALIZED  FINITE  ELEMENT  METHOD 


IVO  BABUSKA,  VICTOR  NISTOR,  AND  NICOLAE  TARFULEA 


Abstract.  We  propose  a  method  for  treating  the  Dirichlet  boundary  condi¬ 
tions  in  the  framework  of  the  Generalized  Finite  Element  Method  (GFEM). 
We  use  approximate  Dirichlet  boundary  conditions  as  in  [12]  and  polyno¬ 
mial  approximations  of  the  boundary.  Our  sequence  of  GFEM-spaces  con¬ 
sidered,  S^,  n  =  1,2,...  is  such  that  S ^  Hq(Q),  and  hence  it  does  not 
conform  to  one  of  the  basic  FEM  conditions.  Let  h M  be  the  typical  size 
of  the  elements  defining  S ^  and  let  u  E  be  the  solution  of  the 

Dirichlet  problem  —A u  =  /  in  Q,  u  =  0  on  dfi,  on  a  smooth,  bounded 
domain  Q.  Assume  that  ll^ll f°r  ^  *-V 

and  | u  —  wj|i/i(f2)  <  C/i™||u|| Hm+1(£i)’  u  ^  Frm_,_1(D)  D  Hq(£ 2),  for  a  suit¬ 
able  uj  E  Then  we  prove  that  we  obtain  quasi-optimal  rates  of  con¬ 

vergence  for  the  sequence  u ^  E  S ^  of  GFEM  approximations  of  u,  that  is, 
|| u  —  Ufj, ||ifi(f2)  <  C h,™ \\ f\\ Hm- i  (Qy  Next,  we  indicate  an  effective  technique 
for  constructing  sequences  of  GFEM-spaces  satisfying  our  conditions  using 
polynomial  approximations  of  the  boundary.  Finally,  we  extend  our  results  to 
the  inhomogeneous  Dirichlet  boundary  value  problem  —A u  =  f  in  D,  u  =  g 
on  dQ. 
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1.  Introduction 

In  the  past  few  years,  meshless  methods  for  the  approximation  of  solutions  of 
partial  differential  equations  have  received  increasing  attention,  especially  in  the 
Engineering  and  Physics  communities.  The  reasons  behind  the  development  of 
such  methods  are  the  difficulties  associated  to  the  mesh  generation,  particularly 
when  the  geometry  of  the  domain  is  complicated.  As  in  the  case  of  the  usual  Finite 
Element  Method,  one  of  the  major  problems  in  the  implementation  of  meshless 
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methods  is  the  enforcement  of  Dirichlet  boundary  conditions.  It  is  the  purpose  of 
this  paper  to  address  the  problem  of  enforcing  Dirichlet  boundary  conditions  in  the 
Generalized  Finite  Element  Method  framework. 

The  classical  Rayleigh-Ritz  methods  for  elliptic  Dirichlet  boundary  value  prob¬ 
lems  assume  that  the  trial  subspace  functions  fulfill  the  boundary  conditions.  Nev¬ 
ertheless,  the  construction  of  such  subspaces  implies  many  difficulties  in  practice 
when  the  boundary  of  the  domain  is  curved.  Therefore,  several  approaches  are 
known  for  dealing  with  the  Dirichlet  boundary  conditions.  One  approach  is  to 
modify  the  variational  principles  by  adding  appropriate  boundary  terms  so  that 
there  will  be  no  need  for  the  trial  subspaces  to  fulfill  any  condition  at  the  bound¬ 
ary.  See  the  works  of  Babuska  [2,  3],  Bramble  and  Nitsche  [13],  and  Bramble  and 
Schatz  [15,  16],  among  others,  for  examples  of  how  this  approach  works  in  prac¬ 
tice.  Another  approach  (used  also  in  this  paper)  is  to  use  subspaces  with  nearly 
zero  boundary  conditions.  This  ideea  was  first  outlined  by  Nitsche  [27]  and  further 
studied  by  Berger,  Scott,  and  Strang  [12]  and  Nitsche  [28]. 

Yet  another  approach  is  the  Isoparametric  Finite  Element  Method  or  IFEM  with 
curved  finite  elements  along  the  boundary.  See  [19]  and  references  therein,  or  [18], 
[20,  22,  23,  24,  29,  30],  among  many  others,  for  more  recent  work  and  applications. 
This  approach  is  typically  used  in  connection  with  a  numerical  quadrature  scheme 
computing  the  coefficients  of  the  resulting  linear  systems.  In  the  applications  of  this 
method,  except  in  special  cases  (such  as  when  Cl  is  a  polyhedral  domain)  the  interior 
Cl I,  of  the  union  of  the  finite  elements  is  not  equal  to  Cl,  although  the  boundary  of 
Cl i,  is  very  close  to  dCl.  That  is,  the  approximate  solution  Uh  is  sought  in  a  subspace 
Vh  C  H^Clh)  and  so,  the  homogeneous  Dirichlet  boundary  condition  u  =  0  on  DCl 
is  “approximated”  by  the  boundary  condition  iq,  =  0  on  dCl^.  In  fact,  Uh  is  the 
solution  of  a  variational  equation  a,h{uh,Vh)  =  ( fh,Vh)h  for  ah  Vh  G  14,  where 
a/l(-,-)  is  a  bilinear  form  which  approximates  the  usual  bilinear  form  defined  over 
H1  (Clh)  x  H1  (Clh) ,  and  fh  €  Vt*  approximates  the  linear  form  Vh  €  14  — ►  /nh  fvhdx, 
where  /  is  an  extension  of  /  to  the  set  Cl ^ . 

Our  approach  has  certain  points  in  common  with  the  isoparametric  method  just 
mentioned  in  the  fact  that  we  are  using  polynomial  approximations  of  the  boundary. 
However,  our  method  does  not  require  non-linear  changes  of  coordinates.  Our 
method  combines  the  approaches  in  the  papers  of  Berger,  Scott,  and  Strang  [12] 
and  Nitsche  [28].  Our  definition  of  the  discrete  solution  is  as  in  [12],  whereas  our 
assumptions  are  closer  to  those  of  [28] .  We  have  tried  to  keep  our  assumptions  at  a 
minimum.  This  is  possible  using  partitions  of  unity,  more  precisely  the  Generalized 
Finite  Element  Method  or  GFEM,  a  method  that  originated  in  the  work  of  Babuska, 
Caloz,  and  Osborn  [8]  and  further  developed  in  [6,  9,  10,  25,  26]. 

Our  construction  is  different  from  the  IFEM  in  that  we  do  not  require  compli¬ 
cated  non-linear  changes  of  coordinates.  Moreover,  our  method  uses  non-conforming 
subspaces  of  functions  and  it  does  not  have  to  deal  with  extensions  over  larger  do¬ 
mains.  It  is  closely  related  to  [11]  which  uses  GFEM  for  elliptic  Neumann  boundary 
value  problems  with  distributional  boundary  data.  The  GFEM  is  a  generalization 
of  the  meshless  methods  which  use  the  idea  of  partition  of  unity.  This  method 
allows  a  great  flexibility  in  constructing  the  trial  spaces,  permits  inclusion  of  a 
priori  knowledge  about  the  differential  equation  in  the  trial  spaces,  and  gives  the 
option  of  constructing  trial  spaces  of  any  desired  regularity.  We  mention  that  the 
GFEM  is  also  known  and  used  under  other  names,  such  as:  the  method  of  “clouds,” 
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the  method  of  “finite  spheres,”  the  “X-finite  element  method,”  and  others.  See  the 
survey  by  Babuska,  Banerjee,  and  Osborn  [5]  for  further  information  and  references. 

Let  us  now  describe  the  main  results  of  this  paper  in  some  detail.  Let  Q  C 
be  a  smooth,  bounded  domain  with  boundary  dLl.  Let  f  £  L2{ fl)  and  u  £  H2( £1) 
be  the  unique  solution  of  the  Dirichlet  problem 

(1)  —A u  =  f  on  O,  u  =  0  on  <912. 

Assume  that  a  sequence  C  H1(  0)  of  test-trial  spaces  is  given  and  define  the 
discrete  solution  u M  £  S p  in  the  usual  way:  B^u^Vff)  =  (/,  v^)  for  all  v M 
(see  Equation  (2)  below).  We  do  not  assume  to  satisfy  exactly  the  Dirichlet 
boundary  conditions,  that  is,  we  do  not  assume  S M  C  Hq(Q). 

Let  us  fix  from  now  on  a  natural  number  m  £  N  =  {1,2,...}  that  will  play,  in 
what  follows,  the  role  of  the  expected  order  of  approximation.  We  shall  make  the 
following  two  basic  assumptions: 

•  Assumption  1,  nearly  zero  boundary  values:  Il'ivll-fU'^an)  <  C^jTlI^vll-HAn) 
for  any  £  S ^  and 

•  Assumption  2,  approximability:  for  any  u  €  HJ+1( £2)  fl  Hq( 12),  0  <  j  <  to, 
there  exists  ui  £  S M  such  that  \u  —  Ui\h^(q)  <  C^IMI/D+qn)- 

These  two  assumptions  are  formulated  in  more  detail  in  Section  2. 

Under  Assumptions  1  and  2,  our  main  approximation  result  (proved  in  Section  2) 
is  the  following 

Theorem  1.1.  Let  S ^  C  H1( f2)  be  a  sequence  of  finite  dimensional  subspaces 
satisfying  Assumptions  1  and  2  for  a  sequence  — >  0  and  1  <  p  <  to.  Then 

the  (unique)  solutions  u  and  u ^  of  Equations  (1)  and  (2),  respectively,  with  f  € 
Hp~1( fl)  satisfy 

llu  -  UA  <  CIi^WuWhp+hq)  <  ChP\\f\\HP-Hn), 
for  constants  independent  of  p  and  f. 

In  Assumptions  1  and  2,  >  0  is  a  sequence  that  goes  to  0.  Intuitively, 

will  play  the  role  of  the  “typical  size”  of  the  elements  in  S^.  However,  in  our 
abstract  setting,  we  are  not  assuming  that  5^  is  constructed  in  any  particular  way. 
Assumptions  1  and  2  are  easy  to  fulfill  with  a  “flat-top”  partition  of  unity  and 
polynomial  local  approximation  spaces.  In  Sections  3  and  4  we  provide  examples 
of  spaces  that  satisfy  Assumptions  1  and  2.  In  Section  5  we  extend  our  results 
to  the  non-homogeneous  Dirichlet  boundary  conditions  case  u  =  g  on  dfl. 

In  this  paper,  we  shall  use  the  convention  that  C  >  0  indicates  a  generic  constant, 
independent  of  p,  which  may  be  different  each  time  when  used,  but  is  independent 
of  the  free  variables  of  the  formulas. 

2.  Approximate  Dirichlet  boundary  conditions 

In  this  section,  we  give  a  proof  of  Theorem  1.1.  We  begin  by  fixing  the  notation 
and  then  we  prove  some  preliminary  results. 

Recall  that  17  C  is  a  smooth,  bounded  domain ,  fixed  throughout  this  paper. 
We  shall  fix  in  what  follows  to.  £  N  =  {1,2, . . .},  which  will  play  the  role  of  the 
order  of  approximation.  We  want  to  approximate  u  with  functions  uIJ;  £  S^,  p.  £  N, 
where  S M  C  U1(H)  is  a  sequence  of  finite  dimensional  subspaces  that  satisfy  the 
Assumption  1  and  2  formulated  next.  Our  first  assumption  is: 
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Assumption  1  (nearly  zero  boundary  values).  There  exists  C  >  0  such 
that 

IKIIffi/2(3n)  <  Ch^\\vJHi{n)  for  any  v M  G  S^. 

So  S/j,  does  not  necessarily  consist  of  functions  satisfying  the  Dirichlet  boundary 
conditions.  Let  |'u|//i(o)  :=  [/n  | Vu^dx]1' 2.  Our  second  assumption  is: 

Assumption  2  (approximability).  There  exists  C  >  0  such  that  for  any 
0  <  j  <  m,  any  u  G  FP+1(f2)  fl  -ffd(fi),  and  any  /x  G  N,  there  exists  ui  G  S M  such 
that 

I u  -  ui\h i(O)  <  C,/i^||w||fl'3+i(n)- 

We  now  proceed  to  the  proof  of  Theorem  1.1.  We  first  need  some  lemmas.  Let 
Q  :  Hl(VL)  — >  H(\  (Q)  be  the  iL1-orthogonal  projection  onto  the  subspace  C 

H1( VL)  of  functions  satisfying  the  Dirichlet  boundary  conditions.  Let  :=  Q(S M). 

Lemma  2.1.  We  have  that  ||rc  —  Q(w)||_f/i(n)  <  Ch™\\w\\Hi(n),  w  G  SM. 

Proof.  Let  us  denote  by  E  :  Hlt2{d fl)  — >  H1{ fl)  a  fixed,  continuous  right  inverse 
of  the  restriction  (or  trace  map)  lL1(fl)  — *  lL1/2(dfl).  That  is,  we  have 

(Ev)\an  =  v  and  ||-Eu||hi(q)  <  C\\v\\Hi/2(dn),  v  G  Hl/2{dTl). 

We  can  chose  such  an  extension  map  E  since  the  restriction  Hl( fl)  — >  lL1/2(9 fl) 
is  continuous  and  surjective.  Consider  now  w  G  Hl{ ft)  arbitrary.  Then  we  '■= 
w  —  E(w\gn)  G  and  hence  the  projection  property  gives 

||w  -  Q(w)||Ff!(n)  <  \\w  -  W£;||ffi(0)  =  ||£,(w|an)||_f/i(f2)  <  11^1/2(0^)  ■ 

If  w  G  5m,  Assumption  2  then  gives  ||w  —  Q(^)||fl’1(n)  <  as 

desired.  □ 

In  what  follows,  we  shall  need  the  following  classical  result  [1,  17]. 

Lemma  2.2.  For  u  G  H1( ft)  there  is  a  constant  C  that  depends  only  on  O  such 
that 

IMIxp(fi)  ^  C[M 2hi(Q)  +  !Mli2(an)] • 

From  this  lemma  we  obtain  that  and  ||fM||_ffi(n)  are  equivalent  norms 

on  S^,  with  equivalence  bounds  independent  of  p: 

Lemma  2.3.  There  exists  C  >  Q  such  that  C~1\v^\h^(q.)  <  ||ty  ll/xqn)  <  Cjv^lx/qo) 
for  all  p  large  enough  and  all  G  S^. 

Proof.  From  Lemma  2.2,  we  have 

\\vv\\ ffqn)  ^  C\]vu\ ffqn)  +  llu/illl2(an)]  ^  C\\vn\ w(ti)  +  ll^llffi/^an)] 

<  Clwlffqn)  +  Ch^Wv^W 2Hi^, 

where  the  last  inequality  is  a  consequence  of  Assumption  1.  Therefore,  for  p  large, 
h n  is  small  enough  and  we  get 

||^||x/i(0)  <  C{  1  -  Ch2™)~1/2 

which  is  enough  to  complete  the  proof.  □ 
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Lemma  2.3  allows  us  now  to  introduce  the  discrete  solution  u M  of  Equation  (1) 
using  the  standard  procedure.  Let  B(v,w)  :=  fn  Vv  ■  \7wdx  be  the  usual  bilinear 
form.  For  p  large,  let  us  define  the  discrete  solution  u ^  G  of  the  Dirichlet 
problem  (1)  by  the  usual  formula 


(2) 


=  /  f(x)vn(x)dx,  for  all  v p  G  S 
Jn 


Let  v  be  the  outer  unit  normal  to  dfl  and  dS  denote  the  surface  measure  on  dfl. 
Similarly,  let  G  S^,  for  p  large,  be  the  solution  of  the  variational  problem 


(3) 


B(wt_l,vl_l)  =  /  d„u(x)vll(x)dS(x),  for  all  ^  G 
Jon 


where  u  is  the  solution  of  Equation  (1).  Note  that  we  need  Lemma  2.3  to  justify 
the  existence  and  uniqueness  of  weak  solutions  u M  and  w^. 

Lemma  2.4.  Let  u  be  the  solution  of  the  Dirichlet  problem  (1)  and  let  and 
be  as  in  Equations  (2)  and  (3).  Then  B^u  —  u^—w^^v^)  =  0  for  all  G  S^;  hence 

I U  -  -  w„\ <  \u  -  u^|ffi(n)  for  all  v M  G 

Proof.  This  follows  from  the  fact  that 


B{u,v^)  =  /  Vu-Vv^dx  =  /  fv^dxP  /  {dvu)v^dS{x)  =  Bfu^  +  Vy), 
J  Cl  J  Cl  J  dCl 


for  all  G  Sfj,. 

We  now  proceed  to  estimate  u M  and 

Lemma  2.5.  For  p  large,  the  solution  uy,  of  the  weak  problem  (3)  satisfies 
(4)  II II h1  (ci)  <  Ch^\\u\\H2^, 

with  C  a  constant  independent  of  p  and  u. 

Proof.  The  relation  follows  from 


□ 


||wjffi(n)  <  c\wn\m(n)  =  CB(wll,wtf)  =  C  duu{x)wtl{x)dS{x) 

J  dCl 

<  C'||9i,'u||i2(9n)||ui/i||i2(an)  <  C\\dvu\\L2(dn)\\wl_l\\Hi/2(dQ) 

<  C,/i™||u||ff2(n)||u;/i||ffi(n). 

Therefore  ||wM||_f/i(n)  <  Ch™\\u\\H2((i),  as  claimed.  □ 

Lemma  2.6.  For  p  large,  the  solution  u M  of  the  weak  problem  (2)  satisfies 

(5)  llw/ill/rqn)  <  C\\u\\H2{n), 

with  C  a  constant  independent  of  p  and  u. 

Proof.  Let  us  first  observe  that  Lemma  2.3  and  then  Lemma  2.4  give 

llW'llifqn)  —  =  CB  {u^ , 

=  C[B{u,U!j)  -  (dvu,u^)dn]  <  C[\B(u,u^)\  +  \(duu,u^)on\] 

<  C||u||ffi(n)||nM||ffi(n)  +  C\\dvu\\L-i(aQ,)\\u^\\Li(dQ,) 

<  Cj|'u||tf2(n)||?qi||ffi(Q)  +  C'/i"l||'u||ff2(n)||u/i||ffi(n). 
Now  it  is  easy  to  see  that  ||uM||ffi(n)  <  Cj|u||ff2(n).  □ 
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Now  we  are  ready  to  prove  Theorem  1.1. 


Proof.  We  shall  assume  p  =  m ,  for  simplicity.  The  proof  in  general  is  exactly  the 
same.  Lemma  2.4  and  the  projection  property,  together  with  Lemma  2.5,  give 

(6)  | u  -  \Hi (q)  <  \u  -  -  w^\Hi(n)  +  |wM| 

<  \u  —  M7|//i(n)  +  Ch™\\u\\H2(u)  <  C'/i™||u||77m+i(n), 

where  for  the  last  line  we  also  used  the  approximation  property  (Assumption  2). 

The  estimate  in  the  ff^-norm  is  obtained  from  Lemma  2.2,  Equation  (6),  As¬ 
sumption  1,  and  Lemma  2.5  as  follows 

llM  —  uu\\h1(Q)  <  C[\u  —  +  llWilU^an)] 

<  Ch™\\u\\Hm+Hn)  +  Ch™\\uJm{n)  <  Ch™\\u\\Hm+Hn). 

The  proof  is  now  complete.  □ 


In  view  of  some  further  applications,  we  now  include  an  error  estimate  in  a 
“negative  order”  Sobolev  norm.  We  let  H~l(Q)  to  be  the  dual  of  Hl(Q)  with  pivot 
L2(Q).  Since  0  is  a  smooth  domain,  can  also  be  described  as  the  closure 

of  C°°  (fl)  in  the  norm 


(7) 


IMIff-qn)  =  sup 
0^0 


|(k,  4>)l*(  o)| 


(Note  that,  in  several  other  papers,  H  l(tt)  denotes  the  dual  of  Hq(LI).) 

Theorem  2.7.  Let  0  <  l,  1  <  p,  and  l+p+1  <  m.  Then,  under  the  assumptions  of 
Theorem  1.1,  the  solutions  u  andu M  of  Equation  (1)  and  Equation  (2),  respectively, 
satisfy 


llu  —  ■*vl|,H-I(n)  <  Chl+p+1\\u\\HP+i(Qj), 
for  a  constant  C  >  0  independent  of  p  and  f  €  Hp~1( fl). 

Proof.  The  proof  of  this  theorem  is  an  adaptation  of  the  usual  Nitsche-Aubin 
trick.  Indeed,  let  us  denote  by  F  £  H1+2(ft)  the  unique  solution  of  the  equation 
—A F  =  <fi,  F  =  0  on  dCl,  for  <fi  £  Hl(Ll)  arbitrary,  non-zero.  Then  there  exists 
a  constant  C  >  O,  independent  of  f>,  such  that  ||F||ff(+2(Q)  <  C\\(j>\\Hi^).  By 
Assumption  2,  there  exists  Fj  e  S ^  such  that  |F  —  Fi |t/i(o)  <  Chlff1\\F\\Hi+2(ny 
An  easy  observation,  which  will  be  used  later,  is  that 


(8) 


l-Fjjifqn)  <  Chlffl\\F\\Hi+2(Q,)  +  |^|/fi(a)  <  C\\(j>\\Hi(ny 
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Using  Lemmas  2.1,  2.5,  and  2.6,  together  with  Assumptions  1-2  and  (8),  we 
obtain 


I U  -  u^\\H-^n)  =  sup - T7— - =  sup 

0/0  0/  o 

|B(u  -  -  Fr)| 


<  sup 
0/ 0 


B{u  -  Up,  F )  +  fgn  Upd^FdS 

fan  u^9uFdS | 


Il0llif'(n) 


s^P  rr~rTj - +  sup 

0/0  Il0ll.f/!(f2)  0/0  ll<Pllff!(n) 


~~  0/0  IMIff'ffi)  0/0 


sup 

0/0 


\uAL2(dn)\\dvF\\L2(dn) 


<  Cl/i^+i+1||M||ffP+i(n)  +  Ch™\\u\\H2{n)  +  Ch™\\u\\H2{n) 

<  Cl/i^+i+1||u||ffp+i(Q), 

by  the  assumption  that  p  +  l  +  1  <  m.  □ 


3.  The  Generalized  Finite  Element  Method 

Our  goal  is  to  construct  a  sequence  Slt ,  /i  =  1,2,...,  of  Generalized  Finite 
Element  spaces  that  satisfy  the  two  assumptions  of  the  previous  section.  To  this 
end,  we  shall  introduce  a  sequence  of  Generalized  Finite  Element  spaces  that  satisfy 
certain  conditions  (Conditions  A(/i;i),  B,  C,  and  D).  In  the  following  sections  we 
shall  prove  that  these  conditions  imply  Assumptions  1  and  2. 

We  begin  by  recalling  a  few  basic  facts  about  the  Generalized  Finite  Element 
Method  [5,  10,  26].  This  method  is  quite  convenient  when  one  needs  test  or  trial 
spaces  with  high  regularity. 

3.1.  Basic  facts.  Let  k  £  Z+.  We  shall  denote  as  usual 

Mw*.°»(n)  :=  max  ||<9Qw||Loc,(n),  ||u||iy*,oo(n)  :=  max  ||dau||Loo(n), 

|  a  | = Ac  |  o:  |  <  A; 

Wk’ 00 (ft)  :=  {u,  IMIm^.oo^)  <  oo},  and  ||  Vw||^,oo(n)  :=  \\dju\\w*-°°(Q)-  In 

particular,  Mwo,oo(n)  =  || ^11  vk°-=° (f2)  =  IMU°°(fi)- 

We  shall  need  the  following  slight  generalization  of  a  definition  from  [10,  26]: 

Definition  3.1.  Let  U  C  R"  be  an  open  set  and  {u>j}jLi  be  an  open  cover  of  fi 
such  that  any  x  €  fi  belongs  to  at  most  k  of  the  sets  u>j.  Also,  let  {<t>j}  be  a  partition 
of  unity  consisting  of  Wm’°°(Q)  functions  and  subordinated  to  the  covering  {u>j} 
(■ i.e supp^>j  C  u>j).  If 

(9)  H^^IUoocn)  <  C'fc/(diamwi)fc,  k  =  \a\  <  m, 

for  any  j  =  1, . . . ,  N,  then  {<fj}  is  called  a  ( k ,  Co,  C i, . . . ,  Cm)  partition  of  unity. 

Assume  also  that  we  are  given  linear  subspaces  'Stj  C  Hm(LUj),  j  =  1,2, ...  ,7V. 
The  spaces  Tj  will  be  called  local  approximation  spaces  and  will  be  used  to  define 
the  space 

N 

s  =  Sgfem  :=  {  Y,  vi  e  c 

3= 1 


(10) 
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which  will  be  called  the  GFEM-space.  The  set  {uij,  <pj,  'Ey}  will  be  called  the  set  of 
data  defining  the  GFEM-space  S.  A  basic  approximation  property  of  the  GFEM- 
spaces  is  the  following  Theorem  from  [10]. 


Theorem  3.2  (Babuska-Melenk).  We  shall  use  the  notations  and  definitions  of 
Definition  3.1  and  after.  Let  {fij}  be  a  (k,Co,C\)  partition  of  unity.  Also,  let 
v j  £  Tj  C  FI1  (uij),  uap  :=  'Yhj4>jvj  €  S,  and  dj  =  diam uij,  the  diameter  of  to j. 
Then 


Xaplll,2(fi) 


<  KCl^\\u- 


Vj  lll2(^)  and 


(11) 


||V(u-uop)||i,2(n)  <  2k  ^  ( 


Cl\\u -v^ 

[dj? 


C'o2I|V(u-^)|||2(^) 


3.2.  Conditions  on  GFEM  data  defining  Recall  that  ui  is  star-shaped  with 
respect  to  w*  C  w  if,  for  every  x  £  ui  and  every  y  £  lo*,  the  segment  with  end 
points  x  and  y  is  completely  contained  in  to.  Also,  recall  that  we  have  fixed  an 
integer  m  that  plays  the  role  of  the  order  of  approximation.  Let  {uij,  f/>7, 
be  a  single,  fixed  data  defining  a  GFEM-space  S,  as  in  the  previous  subsection, 
and  let  E  :=  {uij,  <f>j,  4 !  j,ui*},  where  uij  is  star-shaped  with  respect  to  ui*  C  uij.  We 
shall  need,  in  fact,  to  consider  a  sequence  of  such  data 

(12)  Sp  =  {u%,  <%,  *?,  <*>?}%,  H  £  N, 
defining  GFEM-spaces  Sp 

iV„ 

(13)  *V=  {  E  ’  G  e  ^  }  C  H' m  (fi) , 

i= i 

such  that  there  exist  constants  A,  Cj ,  a,  and  k  and  a  sequence  hp  — >  0,  as  y  — >  oo, 
for  which  EM  satisfies  Conditions  A(/i,J,  B,  C,  and  D  below  for  y  £  N. 

Condition  A (hM).  We  have  that  O  =  U f=iUij  and  for  each  j  =  1,2, . . . ,  Np ,  the 
set  uij  is  open  of  diameter  dj  <  hp  <  1  and  ui C  uij  is  an  open  ball  of  diameter 
>  adj  such  that  uij  is  star-shaped  with  respect  to  ui *G 

Notice  that  we  only  assume  the  open  covering  {uij}  to  be  nondegenerate,  a 
weaker  condition  than  quasi-uniformity  (see  [17],  Section  4.4,  for  definitions  and 
more  information  on  these  notions). 

Condition  B.  The  family  is  a  {n,  Co,  C\, . . . ,  Cm)  partition  of  unity. 

The  following  condition  defines  the  local  approximation  spaces  T)' .  To  formulate 
this  condition,  let  us  choose  Xj  £  uij  (~l  dil.  if  the  intersection  is  not  empty.  We 
can  assume  that  linear  coordinates  have  been  chosen  such  that  Xj  =  0  and  the 
tangent  space  to  dD  at  Xj  is  {xn  =  0}  =  R”_1.  For  /;;t  small,  we  can  assume 
that  uij  fl  dD  is  contained  in  the  graph  of  a  smooth  function  gff  :  Rn_1  — »  R.  If 
x  =  {x\,x2,  ■  ■  ■  ,xn)  £  R",  then  we  shall  denote  x ’  —  {x\,x2,  ■  ■  ■  ,a;„_i)  £  R"_1,  so 
that  x  =  (x',  xn).  Let  qj  :  R"-1  — >  R  be  a  polynomial  of  order  m  such  that 


(14)  |<«)  -<(z')l  <  C{d^)m+1  and 

\X7gj(x')  -  S7qf(x') \  <  C(d^)m  for  all  (x',x„)  £  ui f. 
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This  condition  is  satisfied,  for  instance,  if  dagh( 0)  =  daq^( 0),  for  all  |a|  <  to.  In 
this  case,  the  TO-degree  polynomial  q h  :  Kn_1  — >  K  is  uniquely  defined  by  the  afore 
mentioned  requirement. 

Next,  denote  by  q^  :  R"  — » IR”  the  bijective  map 

(15)  q?(x)  =  q?(x',xn)  =  (x'  ,xn+ q^(x')). 

Let  us  denote  by  Vk  the  space  of  polynomials  of  order  at  most  k  in  n  variables. 

Condition  C.  We  have  T)1  =  Vm  if  D  dLl  =  0  and,  otherwise, 

=  {p  °  (qf  )“\  P  e  Vm,  such  that  p(x',  0)  =  0}, 

where  q h  are  polynomials  satisfying  Equation  (14)  with  a  constant  C  independent 
of  j  and  p. 

An  equivalent  form  of  the  condition  “p  £  Vm,  p{x' ,  0)  =  0”  is  “p  =  xnp\, 
Pi  G  Vm- 1,”  because  any  polynomial  vanishing  on  the  hyperplane  {xn  =  0}  is 
a  multiple  of  xn.  Since  (ql^)~1(x'  ,xn)  =  (x',xn  —  qj(x')),  we  obtain  p(x',xn)  = 

(xn  -  Qj(x'))pi  o  (eft)"1. 

Condition  D.  We  have  (j)h  =  1  on  for  all  j  =  1, . . . ,  for  which  S  0. 

The  constants  Cj,  a,  and  k  will  be  called  structural  constants.  Note  that  we 
must  have  N M  — >  oo  as  p  — >  oo. 

The  above  assumptions  are  slightly  weaker  than  the  ones  introduced  in  [11].  For 
instance,  Condition  C  implies  the  following  propety  (which  is  similar  to  Condition 
C  in  [11]) 

For  any  w  £  Vfd1,  any  0  <  l  <  m  +  1,  and  any  ball  uj*  C  toh  of  diameter  >  adk . 

(!6)  IMIffRup  <  C\\w\\Hl(w,y 

For  further  applications,  we  shall  also  need  a  variant  of  the  spaces  S,,  in  which  no 
boundary  conditions  are  imposed.  Recall  the  functions  q h  used  to  define  the  spaces 
Tj .  Let  =  tEd*  if  LOj  does  not  touch  the  boundary  Oil  and  T)'  =  {po  (qj)^1,  p  G 
Vm}  otherwise,  (the  difference  is  that  we  no  longer  require  p  to  vanish  when  xn  =  0). 
We  then  define 

(17)  ^  ;=  {  E  <%vi  >  vi  e  }  c  Hm  (Si) . 

j= i 

We  shall  also  need  the  following  standard  lemma,  a  proof  of  which,  for  s  £  Z+, 
can  be  found  in  [11].  For  s  >  0  it  is  proved  by  interpolation. 

Lemma  3.3.  Let  ifj  be  measurable  functions  defined  on  an  open  set  W  and  s  >  0. 
Assume  that  there  exists  an  integer  k  such  that  a  point  x  £  W  can  belong  to  no 
more  than  k  of  the  sets  supp (ipj).  Let  f  =  Then  there  exists  a  constant 

C  >  0,  depending  only  on  k,  such  that  ll/ll^q^E)  —  W^iWh^w)- 

Recall  that  dj  denotes  the  diameter  of  (x'J- .  Let  us  observe  that  Condition  A(/im) 
implies  the  following  inverse  inequality. 
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Lemma  3.4.  There  exists  C  >  0,  depending  only  on  a,  such  that 

(18)  Mh*W)  <  C(d^s\\p\\HrK), 

for  all  0  <  r  <  s  <  m,  all  j,  all  p,  and  all  polynomials  p  of  order  m. 

Proof.  The  proof  of  this  lemma  is  inspired  from  the  proof  of  (4.5.3)  Lemma  of  [17]. 
Consider  p  and  1  <  j  <  N/t  arbitrary,  but  fixed  for  the  moment.  Let 


where  x k  is  the  center  of  the  ball  to*h . 

If  p  G  Vm  is  a  polynomial  of  order  m,  then  p  is  defined  by  p(x)  :=  p(djX  +  x^ ) 
for  all  x.  Observe  that  the  set  Vm  :=  {p  :  p  £  Vm}  is  nothing  but  the  set  of  all 
m-degree  polynomials  in  x.  Clearly, 

(19)  \p\h*(q»)  =  {dj)k~n/2\p\HHu>f)’  for  0  <  k  <  m. 

We  first  prove  (18)  for  the  case  r  =  0.  Since  Vm  is  finite  dimensional,  we  have 
by  the  equivalence  of  norms  on  the  unit  ball  B( 0, 1)  that 

(20)  l|p||iffc(B(o,i))  <  C'IIpIU2(b(o,i)),  for  any  0  <  k  <  m, 

where  C  >  0  is  a  constant  that  does  not  depend  on  k,  j,  and  p.  From  Condition 
A (h/f),  we  obtain  that 

(21)  ||p||b2(b(o,i))  <  C\\p\\L2(Q*^ 

where  C  >  0  depends  only  on  the  structural  constant  a.  From  (20)  and  (21),  it  is 
clear  that 

IIpII <  C\\p\\m*f)  Vp  G  Vm, 
where  C  >  0  depends  only  on  a.  Therefore,  (19)  implies 

\p\H<'(w»)(dj)k~n/2  <  C'||p||L2(wp(^)“"/2  for  0  <  k  <  s, 
from  which  we  deduce  that 

\p\h^)  <  C(dj)~k\\p\\L2(u»)  for  0  <  k  <  s. 

Since  dj  <  /i;i  <  1,  we  have 

(22)  Mh-W)  <C(dk)-s\\p\\L,H), 

which  is  just  (18)  for  r  =  0. 

Let  us  now  analyse  the  general  case  0  <  r  <  s  <  m.  For  |a|  =  k,  with  s  —  r  < 
k  <  s,  Dap  =  DdD^p  for  \/3\  =  s  —  r  and  |q|  =  k  +  i —  s.  Therefore, 

H-D“pI|l2(uA‘)  <  ||-D7p||ffa-r(a,p 

<C[$y-'\\D''p\\„{u>)  (by  (22)) 

<  c{qy-'\P\H*+r-.{u»). 

Since 

\p\Hk(ujf)  ;=  X!  \\Dav\\L^), 

\ot\—k 

we  obtain  that 

\p\h^)  <  C(d^)r~s\p\Hk+r-s(^)  for  s  —  r  <  k  <  s. 
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This  implies  that 

(23)  \p\h*(^)  <  C(dj)r~s\\p\\Hr(w»)  for  s-r<k<s. 

From  (22),  we  also  have 

(24)  Ibll *.-rK)  <  C(dj)r~s\\p\\L2^)  <  C{d*y-' \\p\\h^). 

Combining  (23)  and  (24)  gives  (18)  and  this  ends  the  proof  of  the  lemma.  □ 


4.  Properties  of  the  spaces 

In  this  section,  we  establish  some  properties  of  the  GFEM  spaces  Slt ,  p  €  N, 
defined  in  Equation  (13)  using  the  data  E M  =  {ujj,  <j)j,  ,  c satisfying  condi¬ 

tions  A(/im),  B,  C,  and  D  introduced  in  the  previous  section  for  h ^  — »  0.  The  main 
result  is  that  the  sequence  S M  satisfies  Assumptions  1  and  2  of  the  first  section. 

Hereafter,  for  simplicity,  we  will  omit  the  index  p  whenever  its  appearance  is 
implicit. 

Let  us  fix  j  such  that  ZJ jP\dQ  is  not  empty.  Recall  the  functions  fjj ,  q.j  :  Kn_1  — >  K. 
defined  in  the  previous  section.  So,  for  h  small,  ujjCuXl  is  contained  in  {(x' ,  gj(x'))}, 
the  graph  of  the  smooth  function  gj  :  ]R”_1  — >  ffi.  (this  may  require  a  preliminary 
rotation,  which  is  not  included  in  the  notation,  however,  for  the  sake  of  simplicity). 
Let  qj  :  Rn  — »  Kn  be  the  bijective  map  defined  by  Equation  (15).  Similarly,  let 

(25)  gj(x)  =gJ(x',xn)  =  (x1 ,  xn  +  gj(x')). 

Then  gj  maps  K”_1  to  a  surface  containing  ujj  n  dfl.  We  have  gj 1(x)  =  (x',xn  — 
gj{x'))  and  qjl{x)  =  (x',xn  -  qj( x')). 

We  shall  need  the  following  estimate. 

Lemma  4.1.  For  any  polynomial  p  of  order  in,  we  have 

WpogJ1  -P°qJ1\\L^(Uj)  <  Cdy+1\\p\\Hi{u}j)  and 
I \p°~gjl  -poq^Whh^)  <  Cd™\\p\\Hi{Uj), 

where  C  is  a  constant  independent  of  p,  p,  and  j. 


Proof.  By  Taylor’s  expansion  theorem  in  the  xn  variable,  we  have 
pog~1(x',xn)  =  p{ x',xn  -  gj(x'))  =  p{x',xn)  -  gj(x')dnp(x' ,xn)  +  . . . 


+  (-1)' 


k\ 


./^fc 


-d knp{x',xn)  +  . . .  +  (-1)’ 


,9j{x')T 


- d™p(x',xn ) 


and 


P°q3  (x',xn)  =  p{x’,xn  -  qj(x'))  =  p(x',xn)  -  qj(x,)dnp(xf  ,xn)  +  . . . 

+  (-1  )^a^'An)  +  .  .  .  +  (-1  rM^rdrnp{x^Xn). 


k\ 


ml 


Then, 

I pog^ix'iXn)  -P°qJ1(x',xn)\  =  | p{x',xn  -  gj( x'))  - p{x',xn  -  qj{x'))\ 

<  I  g3{x')  -  qj{x')\  ■  \dnp{x',xn)\  +  ...+  | fJj  ^  -  , ,  q°  —  -  •  \d%p(x',xn)\  +  ... 


k\ 

MXT  -  <&(*')' 


\d™p(x' ,xn)\. 
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From  this  and  the  Cauchy-Schwartz  inequality,  we  obtain 


I P°9j  (x',xn)  ~P°q3  1{x',xn)\2  =  | p(x',xn  -  gj{x'))  -p(x',xn  -  qj{x'))Y 

_  12  ,  ,  (9j(x')k  ~qj(x')k)2  k  ,  ^2 


<  m[(gj(x')  -  qj(x'))2dnp(x',xn)2  +  . . .  + 


(fc!)2 

{gAx'T-qjWT? 

(to!)2 


dnP(x',Xny  +  ... 
d™p(x',xn)2}. 


Notice  that  \gj(x')\  <  dj,  for  all  ( x',xn )  £  u>j,  and  because  qj{x')  =  gj(x ')  + 
0(djl+1),  for  all  ( x',xn )  €  u>j,  we  have 

(. 9j(x')k-qj(x')k )2  =  [gk{x,)-(gj{x,)+0{dljl+1))k}2  <  Cdfm+k) ,  for  k  =  lf,. .  ,m, 
which  in  turn  implies  that 

I P°gJ1(x',xn)-poq~1(x',xn)\2  <  Cdfm+1)[dnp(x',xn)2  +  d2d2p(x' ,xn)2  +  ... 

+  dfk-1)d%p(x',xn)2  +  . . .  +  dfm~1) d™p(x' ,xn)2}. 
By  using  the  inverse  inequality  d^WpWfjk^  <  C\\p\\Hi^.),  we  get 


\\po~9j1  -P°~qj1\\l^j)=  [  \pogj\x',xn)-poq7\x',xn)\2dx'dxn 

Jujj 

<  Cdfm+1)  f  [dnp(x',  xn)2  +  d2d2p(xr,  xn)2  +  . . .  +  dfk~1]  dkp(x' ,  xn)2  +  ... 
Jujj 

+  dfm~l)d™p(x',  xn)2]dx'dxn 

<  Ctf^m+r)[\\p\\2Hi^  +  . . .  +  dfk  l)\\p\\2Hk^.-)  +  ■  •  •  +  dfm  ^WpWh^^I 

<Cdfm+1)\\p\\2H1{u}j), 


and  this  completes  the  proof  of  \\pogj  1  -poq^  1|| l*(w})  <  C<f™+x\\p\\Hi(Uj). 

The  proof  of  \\pogJ1  —  po  q3  1\\m(Uj)  <  Cd'jl\\p\\Hi(w:i)  is  reduced  to  the  previous 
inequality  as  follows.  First,  from  the  inverse  inequality  dj ||p||_H-i(w  )  <  C|bl|i2(w -), 
we  obtain 


(26)  \\P°9j  1  ~P°9j  1W^(wj)  <  Cd™\\p\\L2{w.y 
It  is  then  enough  to  show  that 

(27)  \\dkipogJ1)  -dkipoqT1)^^  <  Cd^\\p\\Hi{u)j), 
for  all  k  =  1,2, ...  ,n. 

The  case  k  =  n  is  easier,  so  we  shall  treat  only  the  case  when  l  <  k  <  n  —  1.  A 
Taylor  expansion  with  respect  to  the  ccn-variable  gives 


dk(p  °  9j  1)( x ',  xn)  =  dk(p(x',  xn  -  gj(x'))) 

=  (dkp)(x',xn  -  9j(x'))  -  dkgj(x')(dnp)(x',xn  -  gj(x' )) 

and 


dk(p°qj  1)(x',xn)  =  dk(p(x',xn  -  qj(x'))) 

=  (dkp)(x',xn  -  qj(x'))  -  dkqj(x')(dnp)(x',xn  - 


q3(x')) 
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Equation  (27)  then  follows  from  Equation  (26)  and  from  the  estimates  qj(x')  = 
gj(x')  +  0(djl+1),  dkqj{x')  =  dkgj(x’)  +  0(d J1)  and  \gj(x')\  <  dj  for  (x',xn)  G  u>j 
(see  Equation  (14)  and  Condition  C).  □ 

Remark  4.2.  Let  us  observe  that  Condition  A (h^)  was  used  implicitly  in  the  proof 
of  Lemma  4.1  when  we  used  the  inverse  estimates  d^1\\p\\Hk^w.^  <  C\\p\\ h1^)- 

Remark  4.3.  If  (14)  is  replaced  by  the  more  restrictive  condition  \da(g3  —  qj)  |  < 
C<f^+1~ l“l,  for  all  |a|  <  m+ 1,  then  the  result  of  the  above  lemma  can  be  extended 
as  follows:  For  any  polynomial  p  of  order  m,  we  have 

I \P°9jl  ~P°qJ1\\H‘(uj)  <  CdJ+l~s\\p\\Hi{uJj),  s  =  0, . . . , to  +  1, 
where  C  is  a  constant  independent  of  p,  g,  j,  and  s. 

Corollary  4.4.  Let  p  G  Vm,  then 

Uo{P°9jl  -poqJ^Wm^)  <  Cd™\\p\\Hi{uij). 

If  p  G  Vm  also  vanishes  on  {xn  =  0}  then  we  have 

Il0i(p°971)llffi/2(an)  <  Cd™ |b||ffi(w*). 

Here  C  is  a  constant  independent  of  p,  g,  and  j. 

Proof.  Using  Lemma  4.1  and  Assumption  B,  we  obtain 

UjipogJ1  -poqJ^WHUuj)  <  Uj\\L^(Uj)\\p°9]'1  -poqj'WHUu,,) 

+  ('^i=1\\di(j)j\\Lo.{u!.))\\pogT1  -  pog71||i2(w.) 

<  Cd?\\p\\HHuj)  +  Cd-ldJ+1\\p\\m^)  <  Cdf\\p\\m(uj). 

The  last  part  follows  from  the  first  part  of  this  corollary,  which  we  have  already 
proved,  and  from  the  fact  that  (fj{p°  gj1)  =0  on  Oil.  Indeed, 

Uj(p°qJ1)\\H^(dn)  =  MtipogJ1  --P°<771)ll//i/2(3fi) 

<  C\\(j)j(pogT1  -poqJ1)\\H1{n)  =  CW^^pogT1  -poqJ1)\\H1{u.) 

<CdJ\\p\\m{Wj)<CdJ\\p\\m^) 

The  proof  is  now  complete.  □ 

We  are  ready  now  to  prove  that  Assumption  1  is  satisfied  by  the  sequence  of 
GFEM- spaces  S),  introduced  in  Subsection  3.2. 

Proposition  4.5.  Let  S M  be  the  sequence  of  GFEM-spaces  defined  by  data 
(Equation  (12))  satisfying  conditions  Aih^),  B,  C,  and  D.  Then  the  sequence 
satisfies  Assumption  1. 

Proof.  Let  w3  G^J  and  w  =  4>jWj  G  S^,.  Since  we  are  interested  in  evaluating  w 
at  <912,  we  can  assume  that  only  the  terms  corresponding  to  j  for  which  uj3  n<912  0 

appear  in  the  sum.  Then  Wj  =  Pj  °  qj  ,  for  some  polynomials  pj  G  Vm  vanishing 
on  {xn  =0}.  Hence  Lemma  3.3  and  Corollary  4.4  give 


Miff1/2 (on)  —  II II ff i/2(of2)  —  \\4>j(pj  °  fj  ) II ff i/2(oo) 

3  3 

<  cj2d2r\\P3\\m^)  <  chim  eimhW)  <  chimY,h 


*H1H V 
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By  Condition  D,  ^  ||w,-||^1(wp  =  IMI/ri(Uw;)-  Therefore, 

llu'llffl/2(an)  <  C'/i^m||w||^i(Uu;»)  <  Ch2Jn\\w\\2Hi(ny 
Assumption  1  is  hence  satisfied  by  taking  square  roots.  □ 


Remark  4.6.  Condition  D  is  only  needed  in  the  proof  of  Proposition  4.5.  Although 
one  can  prove  that 

(28)  Hub'll —  CIMlirqn) 

3 

(by  using  norm  equivalence  in  finite  dimensional  spaces),  one  can  not  bypass  Con¬ 
dition  D  because  the  constant  C  in  (28)  depends  on  g .  To  remove  this  dependence, 
one  would  have  to  impose  additional  and/or  different  conditions  on  the  partition 
of  unity. 


The  proof  that  the  sequence  also  satisfies  Assumption  2  is  also  based  on 
the  above  lemma  and  on  the  following  result.  Recall  that  the  local  approximation 
spaces  'Fj  and  4C  were  defined  in  Subsection  3.2. 

Lemma  4.7.  Let  u  £  Hm+1(u>j).  Then  there  exists  a  polynomial  w  £  Tj  such 
that  \\u  -  w\\Hi{ulj)  <  Cdp\\u\\Hm+Huj)  and  || u  -  w\\L^Uj)  <  CdJ+1  ||u||ff™+i(wi) 
for  a  constant  C  independent  of  u,  /i,  and  j.  If  u  =  0  on  uJj  D  dCl,  then  we  can 
chose  w  £  tyj. 

Proof.  We  are  especially  interested  in  the  case  when  u  =  0  on  Uj  D  dCl,  so  we  shall 
deal  with  this  case  in  detail.  The  other  one  is  proved  in  exactly  the  same  way. 

Let  us  consider  v  =  u  o  c/j.  Since  gj  maps  Rn_1  =  {xn  =  0}  to  a  surface 
containing  UJj  C  Oil,  we  obtain  that  v  =  0  on  R™_1.  For  /iM  small  enough,  we  can 
assume  that  g~  (uij)  lies  on  one  side  of  R"_1.  Let  U  be  the  union  of  the  closure  of 
5~1(wj)  and  of  its  symmetric  subset  with  respect  to  Rn_1.  Define  V\  £  Hl(U)  to 
be  the  odd  extension  of  v  (odd  with  respect  to  the  reflection  about  the  subspace 
R”_1).  Let  pi  be  the  projection  of  v\  onto  the  subspace  Vm  of  polynomials  of 
degree  m  on  U .  This  projection  maps  even  functions  to  even  functions  and  odd 
functions  to  odd  functions.  Hence  pi  is  also  odd.  In  particular,  p\  =  0  on  Rn_1. 
We  also  know  from  standard  approximation  results  [17]  that 


IK  _  Pill J?1  (17)  <  Cd™'\\vi\\Hm+i(Uy 


Then 


x— 1 1 


i\u-Pi°9j  1  m(Uj)  <  C\\vi  —  Pi\\m(u)  <  Cd™ 
Let  w  =  pi  o  qf1.  The  lemma  follows  from 


KIKm+l(t/)  <  Cd^WuWHm+i^. 


\\u-w\\Hi{bJj)  <  \\u-p1ogj  1\\m{Uj)  +  \\pi°~9j  1  -P\°9j  1\\m(a,j) 

<  Cdp\\u\\Hm+i(uj)  +  Cd?\\Pl\\Hi(Uj)  <  Cd?\\u\\Hm+iM, 

where  we  have  used  also  Lemma  4.1. 

To  prove  the  relation  ||u  —  io|U2(w3)  <  Cd™+1\\u\\Hm+i(Wjp  we  first  notice  that 
Poincare’s  inequality  gives 

IK  Pi || l2 (u)  <  Cdj\\v i  Pi Wh1  (u)  <  Cd™'+1\\vi\\Hm+nxjy 
The  rest  is  exactly  the  same.  □ 
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We  are  ready  now  to  prove  Assumption  2.  See  [5],  section  6.1,  and  [11]  for  related 
results. 

Proposition  4.8.  The  sequence  of  GFEM  spaces  S M  satisfies  Assumption  2. 

Proof.  We  proceed  as  in  [11],  Theorem  3.2.  Let  u  €  Hm+1(Cl).  If  cZJy  does  not 
intersect  Oil,  we  define  Wj  £  Tj  =  Vm  to  be  the  orthogonal  projection  of  u  onto  Vm 
in  iL1(u’j).  Otherwise,  we  define  Wj  £  using  Lemma  4.7.  Then  let  w  =  JA  4>jwj- 
By  using  Lemma  4.7,  the  definition  of  the  local  approximation  spaces  Tj  (Condition 
C),  and  the  bounds  on  HV^'H^oo  (Condition  B),  we  obtain 

I  u  -  w  |  Hi  (n)  =  |  <t>j  (u  -  wj )  |  (n) 

3 

<  -  Wj\Hi{uj.)  +  ||V^j|z,oo||u  —  Wj\\L2(uj)) 

3 

<  £  (Cd™\\u\\Hm+Hwd)  +  CdJ1dlJl+1\\u\\Hm+i(u,j))  <  CKh™\\u\\Hm+1{u). 

3 

This  completes  the  result.  □ 


5.  Non-homogeneous  boundary  conditions 


In  this  section  we  provide  an  approach  to  the  non-homogeneous  Dirichlet  bound¬ 
ary  conditions.  That  is,  consider  the  boundary  value  problem. 


(29) 


—Au  =  /  on  Cl, 

u  =  g  on  dCl. 


Our  approach  it  to  reduce  again  to  the  case  g  =  0  and  then  to  use  the  results  on 
the  homogeneous  Dirichlet  problem  (1).  In  a  purely  theoretical  framework,  this  is 
achieved  using  an  extension  G  of  g  and  then  solving  the  problem  —  Aw  =  f  +  AG, 
w  =  0  on  dCl.  The  solution  of  (29)  will  then  be  u  =  w  +  G.  This  gives  that 
the  problem  (29)  has  a  unique  solution  u  £  Hp+1(Cl)  for  any  /  €  HP-1  (Cl)  and 
g  £  H1^2+p(dCl)  and  it  satisfies 

IMI/rp+qn)  <  C(\\f\\HP-i^  +  IMIijva+ptan)), 

for  a  constant  C  >  0  that  depends  only  on  Cl  and  p  £  Z+.  (This  result  is  valid  also 
for  p  =  0.) 

In  practice,  however,  we  need  to  slightly  modify  this  approach  since  it  is  not 
practical  to  construct  the  extension  G  (this  is  especially  a  problem  if  g  has  low 
regularity,  that  is,  if  g  is  a  distribution,  for  instance).  We  will  be  looking  therefore 
for  a  sequence  Gk  of  approximate  extensions  of  g. 

The  construction  of  such  a  sequence  of  approximate  extension  as  well  as  the 
analysis  of  the  resulting  method  are  the  main  results  of  this  section.  Other  methods 
for  constructing  Gk  are  certainly  possible.  We  begin  with  an  axiomatic  approach, 
postulating  the  existence  of  the  sequence  Gk,  k  £  N.  Let  1  <  p  <  m.  We  consider 
a  sequence  Gk  satisfying  the  following  assumption.  Recall  that  the  spaces  Sk  A  Sk 
were  defined  in  Equation  (17)  and  are  variants  of  the  spaces  Sk  that  are  not  required 
to  satisfy,  even  approximately,  the  boundary  conditions. 


Assumption  3  (approximate  extensions).  We  assume  that  there  exists  a 
constant  C  >  0  such  that,  for  any  g  £  Hm+1/2(dCl),  there  exists  a  sequence 
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Gk  G  Sk  such  that  ||Gfc|ao  —  g\\H1/2(an)  S  Gh™\\g\\Hm+i/2(Qty  and  ||Gfe||jym+i(Q)  < 

G||(/||fl-™+i/2(an). 

We  now  check  that  it  is  possible  to  choose  Gk  €  A  satisfying  Assumption  3. 
Indeed,  this  is  the  case  since  Assumption  3  is  satisfied  if  the  other  assumptions  and 
conditions  are  satisfied.  We  follow  the  method  in  [4]. 

Proposition  5.1.  There  exist  continuous  linear  maps  Ik  ■  Hm+1(H)  — >  Sk,  such 
that 

(30)  \u  —  /fc(«)|fl-r(Q)  <  Ch™+1  r|M|_F/m.  +  l(Q), 
for  r  =  0  and  r  =  1. 

Proof.  For  u  G  Hm+1(fl)  and  j  fixed,  let  v  =  u  o  g.-.  The  Taylor  polynomial  of 
degree  m  of  v  averaged  over  g~1(ojj)  is  given  by 

(31)  Pj(x)  ■.=  Q™v( x)  =  f  QytVtn(x)^j(y)  dy, 
where 

'  y(a)(y\ 

Qy,v,n{x)  =  v(y)  +  '^2div(y)(xi-yi)  +  . . .+  ^  - — jP-(x-y)a,  a\  =  a  i!...aj, 

i=  1  \a\=m 

is  the  Taylor  polynomial  of  v  at  y  of  degree  m  and  G  C’^{g~l{iOj))  is  a  function 
with  integral  1.  Then,  by  the  Bramble-Hilbert  Lemma,  we  have 

(32)  \v  -  Pj | hs [g~1  (ujj))  <  Gh]?+1~s\v\Hm+i{g-i{w.)),  for  all  0  <  s  <  m+  1. 
Consider  Wj  :=  Pj  o  qj1  G  'F, .  Let  w  :=  <Pjwj-  Then, 


(33)  \u-w\2Hr(n)  <  C^2\<j>j(u-  w)\2Hr{n)  <  C^2\(t>j(u-w)\ 


2 

Hr(vj) 


<cEEi* 

j  i— 0 


j  *= 0 


2  "u-PjO-gr  1'2 


By  changing  variables  and  (32),  we  obtain 


+  \PjO~g-l-Pjo~q-^ 


(34) 


-  Pj  °  9j  =  \vo  ~9j  -  Pj  °  ~9j  <  C\v  ~  Pj 


<  ch2£m+1~r+i'>  wa2 


v\ ~  Ghk  '\UO  9j\Hm.+l(g-l(UJj)) 

<  ni  2(m+l-r+i)|  |2 
-  °'lfc 

Also,  from  Lemma  4.1  and  the  definition  (31)  of  Pj ,  we  have 
(35)  |P,-  057I  -  Pj  o  <  C/l2(-+1-r+i)||p.||^(^) 

<  ChHm+1-r+i)\\u\\Hm+1{uj). 


2(m+l— r+i)  i 
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From  (33),  (34),  (35),  and  Condition  B,  it  follows  that 


(36) 


I u  —  ^ k  2 *[ h 

j  *=o 


— 2irt2(m+l— r+»)  i  i2 


2(m+l— r+i)  i 


U\\  TT 1 


if1  (it 


|j?"*+1(n)> 


for  all  0  <  r  <  m  +  1. 

Define  h-fu)  :=  w.  Clearly  7/c  is  a  linear  map  from  Hm+1(Q)  to  Sk  which  satisfies 
(30).  This  ends  the  proof  of  the  proposition.  □ 


Remark  5.2.  If  we  assume  the  stronger  condition  stated  in  Remark  4.3  on  the  m- 
degree  polynomial  qj,  then  the  conclusion  of  Proposition  5.1  is  valid  for  0  <  r  < 
to  +  1  (the  proof  being  exactly  the  same) . 

From  this  proposition  we  obtain  right  away  the  Assumption  3. 

Proposition  5.3.  For  any  g  G  Hm+1/2(dLl)  there  exists  a  sequence  Gk  G  Sk 
satisfying  Assumption  3. 

Proof.  Let  us  chose  G  G  Hm+1( Q)  extending  g,  ||G||Hm+i(n)  <  C\\g\\Hm+i/2{aa), 
with  C  independent  of  g.  Then  choose  Gk  =  Ik{G),  with  Ik  as  in  Proposition 
5.1.  □ 


Let  Wk  be  the  solution  of 

(37)  — Awk  =  f  +  AGk  on  Q,  Wk  =  0  on  <9f2. 

Also,  let  (wk)fj,  G  Sn  be  the  discrete  solutions  of  this  equation,  namely,  the 
solution  of  the  discrete,  variational  problem 

(38)  B((wk)n,v)  =  (/  +  AGit,i))i2(Q),  v  G 

where  /  G  Hm~1{ f2)  is  the  data  of  Equation  (29). 

The  main  result  of  this  section  is  the  following  theorem. 

Theorem  5.4.  Suppose  Assumptions  1  and  2  and  Conditions  Afk^),  B,  C,  and 
D  are  satisfied.  Let  Uk  '■=  ( Wk)k  +  Gk  G  Sk-  Then  there  exists  a  constant  C  >  0 
such  that  the  solution  u  G  Hm+1{Sl)  of  Equation  (29)  satisfies 

||U  -  Uk  <  ChZ(\\f\\Hm-,{n)  +  ||3||fl-m+l/2(af2)). 

Proof.  We  have  that  Vk  :=  Wk  +  Gk  solves  the  boundary  value  problem 
—A Vk  =  f  on  O,  Vk  =  Gk  on  <912. 

Hence  the  difference  u  —  Vk  solve  the  boundary  value  problem  A (u  —  Vk)  =  0, 
(u  —  Vk)  =  g  —  Gk  on  <912.  From  this  and  Assumption  3  we  obtain 

(39)  ||«  —  UfcHfji(fi)  <  C\\g  -  Gk\\Hl/2(dn)  <  C’^nifl,lljL™+1/2(an)- 
Theorem  1.1  and  Assumption  3,  which  is  satisfied  by  Proposition  5.3,  then  give 

||Dc  —  Wfc||Hi(Q)  =  || Wk  —  (Wfc)fc||ET1(Q)  <  ChZWf  +  AGfcllz/m-qn) 

A  Ch™  (\\f\\Hm-l(n)  +  ||Gfc||/fm+l(fl))  <  C7l™(||/||ffm-l(Q)  +  llplljym+j/a^n)). 


18 


I  BABUSKA,  V.  NISTOR,  AND  N.  TARFULEA 


Hence 

(40)  \\vk  -  Ufellfl-i(n)  =  || wk  -  («Jfc)fc||/fi(n)  <  Ch?(\\f\\Hm-Ha)  +  ||sf||ffm+i/2(9n)). 
The  result  follows  from  Equations  (39)  and  (40).  □ 


6.  Conclusions  and  comments  and  further  work 

We  now  summarize  our  main  results  and  compare  them  with  those  of  Berger, 
Scott,  and  Strang  [12]  and  Nitsche  [28].  We  also  discuss  the  mixed  Dirichlet- 
Neumann  problem  and  suggest  some  further  work. 


6.1.  Conclusions.  Let  p  £  N,  be  the  GFEM  spaces  associated  to  the  “flat- 
top”  data  E^  =  as  in  Subsection  3.2.  In  particular,  they 

satisfy  the  conditions  A(/im),  B ,  C,  and  D,  for  a  fixed  set  of  structural  constants 
A  Cj,  cr,  and  k  and  ht,  — >  0.  (The  parameters  /t;)  are  the  sizes  of  the  patches  u>j.) 
We  know  from  [11]  that  it  is  always  possible  to  find  a  sequence  like  that,  because 
Q  is  smooth. 

Let  ufl  £  S 'p  be  the  discrete  solution  of  the  Dirichlet  problem  (1)  (be.,  —  A u  =  /, 
u  =  0  on  the  boundary  of  the  smooth,  bounded  domain  fi).  That  is,  u M  is  given  by 
Equation  (2).  Then  we  have 

Theorem  6.1.  The  sequence  S f  satisfies  the  Assumptions  1  and  2  and  hence 
llu  ~  un llffi(n)  <  Ch^\\u\\Hm+i^, 

where  C  >  0  is  a  constant  that  is  independent  of  p  and  f  £  iLm_1(H). 


We  therefore  obtain  quasi-optimal  rates  of  convergence  for  our  approximate  so¬ 
lutions  £  S/j, ,  in  the  sense  that  the  error  estimate  has  the  same  order  as  the  best 
approximation  in  the  spaces  S^. 

The  definition  of  the  discrete  solution  u ^  is  as  in  [12].  In  that  paper,  Berger, 
Scott,  and  Strang  obtain  for  to  >  3/2  the  estimate 

(41) 


\u-Uu.\\h1(CI)  <  C 


W- £11^1(0)  +  llM 


Uu\\  jjm(n)  sup 


||u||ff3/2-m(Q) 

IkllffHfi) 


Our  result  is  thus  an  extension  of  the  result  of  [12]  to  the  case  to.  =  1,  to  which  the 
methods  of  that  paper  do  not  seem  to  apply.  The  case  to.  >  3/2  seems  not  to  be 
enough  to  provide  optimal  rates  of  convergence  directly. 

The  definition  of  the  discrete  solution  Uh  €  Sk  in  Nitsche’s  paper  [28]  is  such 
that  (A Uh  +  f,Vh )  =  0  for  all  Vh  £  Sh,  so  it  is  different  from  ours.  In  addition  to 
our  approximability  assumption  (Assumption  1),  in  [28],  Nitsche  also  requires  an 
approximation  property  at  the  boundary  and  ||u||^i(9n)  <  C,^_7||mI|l/1(0)  instead 
of  our  Assumption  2.  These  assumptions,  slightly  stronger  than  ours,  also  lead  to 
optimal  rates  of  convergence. 

As  in  the  isoparametric  methods,  our  method  uses  polynomial  approximations 
of  the  boundary.  However,  we  do  not  have  to  use  non-polynomial  approximation 
functions  of  non-polynomial  changes  of  coordinates. 
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6.2.  Mixed  boundary  value  problems.  Let  us  comment  now  a  little  on  the 
mixed  boundary  value  problem: 

du 

(42)  —A u  =  f  in  12,  u  =  0  on  9dI2,  and  —  =  0  on  9;vf2, 

on 

where  9f2  is  the  disjoint  union  of  the  closed  subsets  9/jf2  and  d^Tl.  In  particular, 
both  9of2  and  9^12  will  be  closed  surfaces.  In  case  9d12  is  empty  (the  Neumann 
problem),  we  require  our  approximation  spaces  S M  to  be  such  that  jn  Vfj,(x)dx  =  0 
for  all  Vfj,  G  S^. 

Then  our  results  remain  the  same  if  we  consider  the  following  variant  of  our 
Assumptions  1  and  2,  in  which  90  was  replaced  with  9^0. 

•  Assumption  lm:  ||tvl|ffi/2(aDn)  <  Ch^\\v^\\Hi(a)  for  any  v p  G  S^,  and 

•  Assumption  2m:  for  any  u  G  LP+1(f2),  0  <  j  <  m,  u  =  0  on  9dO,  there 
exists  ui  G  Sjx  such  that  |u  —  it/li/qn)  <  C^ullull//j+1(n)- 

The  definition  of  discrete  solution  is  the  same  as  before:  u M  G  Stl  is  such  that 

(43)  (V(u  -  Up),  VuM)  =  0,  for  all  G 

An  inspection  of  our  arguments  used  in  the  proof  of  Proposition  4.8  shows  that, 
if  only  Conditions  A(/im)  and  B  are  satisfied,  then  we  can  take  S M  =  {up  G 
Sjj, ,  Jn  uM(x)dx  =  0}.  The  analog  of  Theorem  1.1  then  follows  as  before. 

Theorem  6.2.  The  (unique)  solutions  u  and  of  Equations  (42)  and  (43),  re¬ 
spectively,  with  f  G  Hp~1(  12)  satisfy 

llu  _  uv llffi(n)  <  C'^llulliLp+1(0)  <  Chp\\f\\Hi>-i(n), 

for  constants  independent  of  p  and  f. 

Let  us  observe  that  for  the  mixed  Dirichlet-Neumann  boundary  conditions,  Con¬ 
ditions  C  and  D  can  be  weakend  by  replacing  912  with  9/5 12.  Thus,  in  the  case  of 
pure  Neumann  boundary  condition  case,  Conditions  C  and  D  are  not  required. 
However,  in  the  analysis  of  elliptic  boundary-value  problems  with  mixed  Dirichlet 
and  Neumann  boundary  conditions,  the  Conditions  C  and  D  have  to  be  considered 
only  for  the  indices  j  corresponding  to  the  patches  c 0j  that  touch  the  portion  of  the 
boundary  on  which  the  Dirichlet  boundary  conditions  are  imposed.  See  [4,  5,  11] 
and  references  therein  for  more  information  on  meshless  methods  for  Neumann 
boundary  conditions. 

6.3.  Further  problems.  In  spite  of  all  these  differences  in  assumptions  and  defi¬ 
nitions  between  [12,  28]  and  our  paper,  the  main  issue  seems  to  be  providing  simple 
examples  of  spaces  satisfying  the  various  assumptions  used  in  these  papers.  For 
instance,  it  would  be  interesting  to  provide  other  examples  of  spaces  satisfy¬ 
ing  Assumptions  1  and  2.  It  would  also  be  interesting  to  see  if  a  modification  of 
the  uniform  partition  of  unity  can  give,  by  restriction,  spaces  S ^  satisfying  these 
Assumptions.  Finally,  it  would  be  important  to  integrate  our  results  with  the  is¬ 
sues  arising  from  numerical  integration  and  to  provide  explicit  numerical  examples 
testing  our  results. 

Some  numerical  tests  and  related  theoretical  results  can  be  found,  for  example, 
in  [7]  and  [21]. 
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